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Abstract 

The dissociation energy (Do) of CN is determined to be 7.65±0.06 eV. Th» 
corresponds to AH,(CN)=105.3±1.S kcal/mo e, in excellent agreement with E»- 
gleman and Rouse (1975), but considerably larger than the recent value de uce 
from shock-tube studies by Colket (1984). The above result is obtained not only 
from extensive ul initio MRC1 calculations using a very large Gaussian basis set, 
but also from extrapolation of the directly computed value by comparison of com- 
puted and experimental results for NO, C, aid NV As an additional cal, brat, on 
of our methods, we compute the D„ value for CN from the corresponding value for 
CN- using the experimental electron affinity data. The lifetime of the e = 0 level 
of the violet (B 2 £ + - X 2 £ + ) system has been computed to be 62.4 ns, in i goo 
agreement with both experiment and previovs calculations. Lifetimes for the red 
U 2 n - X 2 £ + ) system decrease with increasing v', which is consistent both wit 
recent experiment and calculations. While our computed lifetimes are 
longer than those obtained from experiment, they are shorter than those e uce 
from an analysis of the solar spectrum. Hov/ever, our Do and /oo are consis en 
with Lambert’s model for the solar spectrum 
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I Introduction 

The CN radical is present in flames, com.ts, stellar atmospheres and inter- 

stellar space. It can be observed by either the red (rh 2 n - X E ) « ™ e 
- X 2 £+) emission systems. In principle, these two CN band systems 

should be an excellent probe of the chemistry u ider these condemns, where erect 
measurements are difficult or impossible. If the dissocial, on energy (Do) and th 
oscillator strength (/oo) - CN (and the temperature) were known the , 0 served 
emission intensity in a stellar atmosphere could be converted into 
dances. Considering the importance of CN as a probe spaces, rt is ““ 
that, while the oscillator strength of the violet system is well known, **■>*•» 
suits exist for both the oscillator strength of the red system and the 8 ™" 
dissociation energy. Since the red system is especially important as a n troge 
abundance indicator for red giants, and contr butes to the total opacty of a red 
giant’s atmosphere (Sneden and Lambert 1982), it would be of enormous valu 
the understanding and modeling of stellar atmospheres to have accurate 
both the ground state dissociation energy and he oscillator streng o 
tern. Such a goal is difficult to accomplish by experimental means - the absence 
suitable predissociations complicates the dete urination of Do, ^ 
surements for the red system are substantial!; affected by (non-ra ra ive 
intersystem crossing. This latter effect is qmte large, leading to the observe . - 
crease in emission with increase in pressure up to 0.02 torr (Dune, Erman an 
Larsson 1978). Indirect methods for determining the rl’Il state lifetime epen on 
a knowledge of Do, so the uncertainty in D 0 further hampers an understandmg o 

the red system. 

From the standpoint of 06 initio quantum chemistry, CN is no more difficult 
than other first-row diatomic molecules wit;, a similar numberofelectrons. 
have recently presented (Langhoff, Bausch, cher, and Taylor 1987; Bau.chhcll^ 
and Langhoff 1987a) very accurate spectros :opic constants or a num 
molecules; advances in the methodology of electronic structure calculations have 
made it possible to compute, for example, the D„ values for C„ N 2 Or and , 
with an error of less than 0.18 eV in the worst case. By using such results o pro 
a calibration of techniques for extrapolating . he best computed value to expenmen , 
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the calculated CN D 0 values can be corrected for the remaining errors arising from 
truncation of both the one- and n— particle besis set expansions. 

In addition to using the C 2 and N 2 moleci les for calibration, we also compute 
Do values for NO and CN - . The NO value is i sed to verify that the extrapolation 
procedure works well for a heteronuclear diaton ic, while an independent calibration 
is obtained by using the computed Dq of CN - h conjunction with the experimental 
electron affinities (EAs) of CN and C (Hoto;> and Lineberger 1975; Berkowitz, 
Chupka and Walter 1968), via the relation 

D 0 (CN)+EA(CN)=Do(CN- )+EA(C). (1) 

Ideally, the use of equation (l) with the compvted D 0 for CN - will yield the same 
Do for CN as that computed (or extrapolated ; directly. An advantage of equation 
(1), however, is that CN - dissociates to C - ( 4 S) and N ( 4 S), which avoids the 
need to describe 2s 2 -» 2p 2 near-degeneracy ejects in the C{ 3 P) state. 

Recent calculations of spectroscopic intensities also show excellent agreement 
with experiment (see, for example, Bauschlicher and Langhoff 1987b), and such cal- 
culations can readily be performed for CN. The availability of good experimental 
data for the violet system (Duric, Erman and Larsson 1978; Jackson 1974) provides 
an additional calibration of the theoretical results for the red system where there is 
a wide variation in the experimental data. It is pertinent to note that computed life- 
times for several states of Nj (Langhoff, Bauschlicher and Partridge 1987, Langhoff 
and Bauschlicher 1987), which is isoelectronic with CN and has similar A - X and 
B - X separations, are in excellent agreemen with experimental results (Holland 
and Maier 1972; Peterson and Moseley 1973). 

In view of the problems that the absence of reliable CN data generates for as- 
tronomical modeling, and given the accuracy cf present large-scale ab imtio calcula- 
tions, we have determined the Do for the CN g ound state and radiative lifetimes for 
the A 2 II and £ 2 E+ states. These results are compared with current experimental 
estimates and with previous theoretical calcu ations. 

II. Methods 

The one-particle basis sets employed in th is study are constructed using general 
contractions based on atomic natural orbitals ANOs) as described in previous work 
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(Almlof and Taylor 1987). We use two different contractions, [4s 3 p 2d If] and 
[5s 4p 3 d 2/ 1 y], of the (13s 8p 6 d Af 2g) primitive Gaussian set. The s and p 
exponents are taken from van Duijneveldt (19 "l) and the polarization functions 
are given in Bauschlicher, Langhoff and Taylor >987) and Langhoff, Bauschlicher, 
and Taylor (1987) for C and N, respectively. The difference in the computed D 0 
values obtained with these two contracted bas s sets provides a measure of basis 
set incompleteness. For CN~, the basis sets art supplemented with a diffuse p set 
on C and N (o p (C) =0.026 and a p (N) =0.038) m describe the more diffuse charge 
density of the negative ion. Finally, an extended [5s 5p Ad 3/ 1 g] ANO set is used 
as an additional test of basis set saturation. Pure spherical harmonic functions are 
used in all calculations. For the calculation of the lifetimes of the red and violet 
band systems, potential energy curves and transition moments are computed using 
the smaller [4s 3p 2d If] ANO basis set, as calibration calculations show that the 
lifetimes computed using this basis are essential y the same as those using the larger 

[5s 4p 3 d 2/ lg] set. 

In this study we use multireference configuration-interaction (MRCI) wave 
functions that are based on a complete-active-* pace self-consistent field (CASSCF) 
zeroth-order reference function. We also use the modified coupled-pair functional 
(MCPF) method (Chong and Langhoff 1986; Ahlrichs, Scharf, and Erhardt 1985), 
based on SCF orbitals, in investigating the oie-particle basis effects on D 0 . For 
CN _ , NO and N 2 , the 2p electrons and orbitals are treated as active in the CASSCF 
calculations. For CN and C 2 the 2s and 2p o bitals and electrons are active: this 
is necessary to account for the 2s-2p hybridize tion in the molecules and 2s 2 -» 2p 2 
excitation in the S P state of C atom. In the study of the A 2 I1 and £ 2 £ + state life- 
times, a single set of CASSCF orbitals is used, obtained from a calculation in which 
the average energy of the X, A and B states is optimized with equal weighting for 
the three states. In all MRCI (and MCPF) calculations, the 2s and 2 p electrons 
are correlated and the two la-like core orbitAs are constrained to be doubly occu- 
pied. Since the configuration expansions including all single and double excitations 
from all configurations in the CASSCF are very lengthy, in many calculations only 
selected references are included. For CN, the MRCI treatment with all CASSCF 
configurations used as references correspond* to a second-order Cl (SOCI) calcu- 
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1 at ion. The selection of references is calibrated igainst SOCI calculations in the 
[4s 3 p 2d If) basis set. In the selection of references, an occupation is included in 
the reference list if the absolute value of the coefficient of any one of its component 
spin couplings is above a designated threshold. We construct reference lists either 
by using only geometries close to r e and at long range, or by merging the lists 
obtained at various points over the full potential curve. The former approach is 
used in the calculation of D 0 , while the latter is used for the determination of the 
radiative lifetimes. In order to estimate the imp ortance of higher excitations, we 
use the multireference analog of the Davidson c< rrection (Langhoff and Davidson 
1974; Blomberg and Siegbahn 1983) denoted +Q and defined as AE(1 - Y,r c r)i 
where A E is the MRCI correlation energy and Cr is the coefficient of reference 

configuration R in the MRCI expansion. 

The computation of radiative transition pro babilities involves the expectation 
value of the electronic transition dipole operator between the X 2 E + , A 2 U and 
B 2 L + states. For the perpendicular ( A -* X) ransition, the transition moment 
operator, D(r), is defined as 

D(r) — < X 2 E + 1^— ^ — |j4 2 II > (2) 

v2 

Note that this definition of the matrix element is equivalent to the cartesian form 
< X 2 Z+\x\A 2 n z > (Whiting tt al ., 1980). For the parallel (B -» X) transition, 
the moment is the usual 

D(r) = < X 2 E^\z\B‘E + > . (3) 

To compare with experimental observables such ;is lifetimes, the transition moment 
operator is integrated numerically over the nuch ar coordinate r, 

rOO 

D v >y = / Vv(r)£>(r)^v" ( r ) dr > 

Jo 

where rp v ' and rp v » are the upper and lower state vibrational wave functions. Note 
that rotationless potentials are used throughout, but the influence of vibration- 
rotation interactions on the Einstein coefficients and lifetimes is rather small. The 
vibrational eigenfunctions are determined using Klein-Dunham potentials and the 
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theoretical transition moments. The Klein-Dunham potentials for the X 2 E+ and 
A 2 II states were based on the Dunham coefficients given by Kotlar et al. (1980), 
while the B 2 £ + potential was based on the constants given by Huber and Herzberg 
(1979). Since the theoretical potentials are in excellent agreement with these Klem- 
Dunham potentials, the Einstein coefficients and radiative lifetimes would not be 
substantially different if the vibrational analysis vas done using entirely theoretical 
input. The MRC1 T t values (based on a 0.2 a< grid of points) for the A 2 II and 
B 2 £ + states (with experimental values in parentheses) are 9560(9245) cm and 
26336(25752) cm -1 . This indicates that the CASSCF/MRCI procedure provides 
a balanced treatment of the three states. As no ed above, the transition moments 
were computed using the smaller [4s 3p 2d If] basis and a 0.05 selection threshold 
on reference CSFs. The vibrational eigenfunctions in equation (3) were obtained 
from solutions of the radial Schrodinger equation using the numerical integration 
method of Cooley (1961). 

Once the dipole matrix element is known, tl e Einstein A coefficient (s ) given 

as 

A v > v » — 2.026 x 10 6 <7 f v'v" I 

can be computed, where P v > v » is the transition energy in cm -1 between levels v' 
and v", and g is a statistical weighting factor, equal to unity for II -► £ and £ -> £ 
and transitions, but two for £ -» II transitions. The radiative lifetime (v) of 
a vibrational level depends on the sum of the .ransition probabilities to all lower 
vibrational levels in all lower electronic states. If rotational effects are neglected, 
the lifetime can be written as 



All calculations were performed using the MOLECULE-SWEDEN codes 
(Almof; Siegbahn et al.) on the NASA Ames CP AY X-MP/48 and the NAS Facility 

CRAY 2. 

III. Results and Discussion 
A. Calculation of Do 

One question raised in a previous study of the dissociation energy of CN by 
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QWbahn and Aeren (1983) is whether the N 2s correlation must be in 
eluded As can be seen in Table 1, excluding th, N 2s orb.ta, Iron, the CASSCF 
active space decreases the Do by 0.34 eV, most of which (0.21 eV) carnes over 
MRC, level Thus a nine-electron treatment is . ecessary for q uant,ta„ve results. 
While it may seem contradictory that including N 2s correlation ° 

CN but decreases Do for Nr (Langholf, Bauschlicher and Taylor '^7), the. wo 
cas K are not, in fact, parallel. In CN the third c orbital is the N 2s orb.tal on y 
large r. at r, it is significantly delocalised. Thu, as the bond length var.es m. 
different fractions of N 2s and C 2s correlation are recovered, an it .s no p 
to relate this to Nr, where eliminating 2s correlation has a much more cons.ste 

effect across the potential curve. 

Since the SOC1 expansions for CN in the 4p 3d 2/ . 9 ] basis ^ 

v h ;ii; An pcFk we have considered the effect ol reter 

electrons correlated are over 7 rn.ll.on CSFs w ^ ^ 

ence selection on the the spectroscope constan-s (Do, r. and .) ^ 

that, although the efTect on r t and w, of reference se ect.on .s sm , 

in Do is considerably larger. The MRC1 and MRCI+Q results at 0.01 referen 

selection threshold in the |4s 3p 2d 1/) basis are compared with 

SOCI and SOC1+Q results in Table 1. Since th.se differences are w.th.n acceptabl 

, ■ , . , ... expansion by about a factor of two, we have 

limits, and selection reduces the Cl expans, on ny 

employed the MRCI(O.Ol) procedure for calculations invo v.ng e | 

basis. However, at 0.025 and 0.05 selection thresholds the C. truncat.on errors 

too large to reliably state error bars to within 0.1 eV. 

Our calculated Do values are reported in Table 2, the C 2 

from Bauschlicher and Langhoff (1987a). Fo. CN - , Nr and NO, al of te- 

figurations in the CASSCF are references in he MRC1 ealeu at.ons (no 

selection). For many of the systems, we also .eport Do values at he MCPF ■ 

As MCPF is a single-reference method, it al ows us to consider the effect 
As M^rr 6 . There are two sources of 

ther basis set saturation at a relat.vely monest cost. There 

error in our calculations, incompleteness of tl e one-particle space (bas,s se, error) 
n-partide space (error in the ccrrelation — , 

(Bauschlicher and Langhoff 1987c) involving full Cl compansons (that .s, ehmina 
tion of the n-partide space incompleteness) indicate that the largest error m Do 
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comes from basis set incompleteness. It is difficult to conceive of a source of error 
in the n-particle space treatment which is pres, nt for CN and CN - and absent for 
the systems studied previously. At the MRCI level, it is therefore vital to reduce 
the one-particle error as far as possible: the residual one-particle and n-particle 
incompleteness can then be accounted for with little uncertainty by extrapolation 
approaches based on the well-knowm D 0 value of N 2 . We consider first the effect 
of expanding the basis set from [4s 3 p 2d If] to [5s 4p 3 d 2/ Iff]. (As noted above 
the C and N basis sets include an extra p set for CN - ). The effect of this basis 
set improvement is very similar for CN, CN - and NO; it is only slightly smaller for 
C 2 , and is somewhat larger for N 2 . For CN and NO, the MCPF calculations show 
the same increase in Do with basis set improviment as do the MRCI calculations. 
For N 2 the effect of basis set improvement is H rger at the MRCI level than at the 
MCPF level, suggesting a coupling of one-parti :le and n-particle errors in N 2 that 
is either not present in the other systems or is b itter reproduced at the MCPF level. 
These results indicate that the remaining errors in CN (or CN - ) should be smaller 
than those in N 2 . Further, for CN the expansion of the basis set to [5s 5p 4d 3/ Iff] 
results in only a small increase in D 0 at the MCPF level, supporting the hypothesis 
that any remaining basis set incompleteness in CN is small. From these observa- 
tions we expect the error in the computed Do for CN will be slightly larger than 
that for C 2 , but smaller than for N 2 . The effect of the +Q correction is small for all 
molecules except NO. This in turn suggests thet the n-particle error in CN will be 
smaller than in NO. The overall similarity of hese systems suggests that NO, C 2 
and N 2 are reasonable choices for the calibrat on of calculations on CN and CN . 
This is in accordance with chemical intuition: CN, N 2 and CN - have a bond order 
of 3, NO has bond order of 2.5 and the bond o? der in C 2 is probably best described 
as between 2 and 3. The energy required for the sp hybridization of C, however, 
results in smaller D 0 values for CN and C 2 thin is indicated by their bond orders. 
This is an important reason for including the C 2 molecule in the calibration study. 

A particular advantage of ANO basis sets is the systematic way in which the 
computed results improve as further ANOs ere added. We may expect that the 
increase in D 0 observed on expanding the basis from [4s 3p 2d If) to [5s 4p 3 d 2/ Iff] 
reflects inclusion of the same effects that wou d be accounted for if the latter were 
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extended to completeness: no new effects are expected to arise. Consequently, if we 
assume that the remaining one- and n-particle errors scale in the same way, we can 
extrapolate to the “experimental” D 0 for CN i sing the well-known experimental 

Do value for N 2 

Po([5s 4 p 3d 2/ lg]) - Pq([4s 3 p 2d 1/;)(CN) 

£> 0 ([55 4 p 3 d 2/ lg)) - I>o(! 4s 3 p 2d If )(N 2 ) 

Do(expt) - Po (l5s 4p 3d 2/ lg])(CN) 

" P 0 (expt) -D 0 (l5s 4p 3d 2/ 1j])(N 2 ) ’ 

If instead, the remaining error is assumed to be proportional to Do as suggested by 
Larsson, Siegbahn and Agren (1983), we obtain the relation 

P 0 ([5s 4p 3d 2/ lg]) _ P 0 (jSs 4p 3d 2/ lg]) 

Po(expt) I>o(expt) 

We refer to the extrapolation formulae (7) and (S) as methods I and II, respectively. 

One test of the accuracy of these extrapolation methods is to apply them to 
the C 2 and NO molecules. The comparison for C 2 is complicated somewhat by the 
remaining uncertainty in the D 0 value of this molecule. Huber and Herzberg (1979) 
report a value of 6.21 eV, which is the average of two thermochemical values. The 
value that we adopt for comparison here is the most recent thermochemical de- 
termination of 6.17±0.11 eV by Kordis and Gingerich (1973). Lower values of Do 
based on spectroscopic methods have also bee r reported (see Lambert (1984) for 
a good review of experimental work). The mest accurate of these is the value of 
6.11±0.04 eV determined by Messerle and Krauss (1967) on the basis of studies of 
the rotational predissociation in the c l U g state. However, since our best directly 
computed value is 6.11 eV and remaining basis set saturation is expected to increase 
this value, we believe that the spectroscopic vdues are slightly too low. Using the 
value of 6.17 eV as our comparison standard, the errors in the C 2 predictions are 
MRCI(MRCI+Q) +0.03 (-0.02) for method I and +0.05(0.00) eV for method II. 
Hence, the +Q values are in perfect agreement with the value of 6.17 eV, while 
those without the +Q correction favor a slightly larger value of Do for C 2 . 

For the NO molecule an accurate Do value is available for comparison see 
discussion in Huber and Herzberg (1979). For NO the errors are -0.12(-0.0l) eV and 
-O.ll(-O.Ol) eV for Methods I and II, respectively. The largest errors are therefore 
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for NO at the MRCI level. However, the Davidson correction effect on Do is quite 
large for NO, which is different from the situation in N 2 . There is thus some 
difference in correlation effects between NO and the calibration standard (N 2 ), and 
it is not surprising that the extrapolation procures fail to account for an effect 
that is unique to NO. The application of the extrapolation methods to CN yields 
D 0 values of 7.64(7.61) eV and 7.67(7.65) eV. Applying the same extrapolations to 
CN“ and computing the D 0 of CN using the experimental EAs and equation (l) 
yields 7.61(7.62) eV and 7.68(7.68) eV. The average of these extrapolated values is 
7.65 eV, leading to a recommended value of 7.65±0.06 eV for CN. The error bars 
correspond to estimated confidence limits of 80%, or approximately 1.5 times the 
average of the above uncertainties. 

By increasing the error bars further, we can increase the confidence limits 
to 99% . First, in order to obtain an upper bound to Do, we assume that the 
maximum error in the computed Do for CN is 1.25 times that for N 2 : this should 
be a conservative upper bound, because there b a smaller change in D 0 with basis 
set improvement and a smaller coupling of basis set and higher excitations for 
CN than for N 2 . In estimating a lower bound, it should be noted that since the 
MRCI treatment used has virtually no Cl superposition error: D 0 could then be 
too large only if the +Q correction substantial!} underestimated the effect of higher 
excitations (note that the +Q correction reduces Do for CN). This is highly unlikely 
as experience has shown that the correct answer for treatments of this type lies 
between the MRCI and MRCI+Q levels (Bauschlicher and Langhoff 1987c), and 
hence the largest reduction from such a source would be the 0.05 eV observed at the 
[4s 3 p 2d If] SOCI level. It is inconceivable tha the remaining basis set errors are as 
small as 0.05 eV, given that remaining basis set error in Do for CH, a simpler system, 
is 0.03 eV (Bauschlicher and Langhoff 1987d), and thus basis set incompleteness 
would at the least cancel any reduction in Do h orn higher excitations. Furthermore, 
the directly computed value based on D 0 for CN“ is 7.50 eV, and in this case 
the +Q correction increases D 0 . Accounting for a minimum remaining basis set 
incompleteness error of 0.05 eV in this value suggests that 7.55 eV would be an 
absolute lower limit to D 0 for CN. This leads us to a very conservative estimate 
(99% confidence limits) of Do 7.65±0.10 eV. 
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A summary of the Do determinations pric r to 1973 is given by Arnold and 
Nicholls (1973); the range of values is 7.50-8.40 eV. These authors noted that a 
range of 7.8-8.0 eV was consistent with their va ue (7.89±0.13 eV) and most of the 
other data available at that time. This range is outside our conservative estimate. 
However, the determination of D 0 by Englemar and Rouse (1975) (7.65±0.05 eV) 
is significantly smaller than that of Arnold and Nicholls (1973), and in excellent 
agreement with our value. While the Engleman and Rouse (1975) determination 
is an indirect result based on the Arnold and Nicholls oscillator strength of the 
violet system and some assumptions about the kinetics of the system, we conclude 
that it is the most reliable experimental deter mination in light of our calculated 
values. Lambert et al. (1986, see also Lambert et of. 1984) obtained an estimate 
of £>o=7.55±0.05 eV from solar red system Hues with the C and N abundances 
obtained from either the CH A - X and NH A - X or CH A - X and the permitted 
N I lines. Their estimate is based on the foo for the red system calculated by Larsson, 
Siegbahn, and Agren (1983). If the smaller f- value obtained in this work is adopted, 
Lambert et al .' s (1986) D 0 value is increased to A>=7.69±0.05 eV (Lambert 1987). 
Our estimates for D 0 and the red system lifetimes are thus consistent with the solar 
data and the model used by Lambert. 

As noted at the beginning of this section, Larsson, Siegbahn and Agren (1983) 
have reported a Do value for CN using a C A SSCF /externally contracted Cl ap- 
proach. Their be6t basis set is somewhat inferior to the smaller [4s 3p 2d 1/] set 
used in the present work, and they select reference occupations using a 0.05 thresh- 
old. Their Do of 7.19 eV is slightly smaller than our SOCI value of 7.32 eV in the 
[4s 3 p 2d 1/] set. Using method D for extrapolation and a comparable calculation for 
N 2 led to an extrapolated value of 7.52 eV, which is smaller than even our directly 
computed result. It is very likely that their basis set is too small and their reference 
selection threshold not sufficiently strict to allow an accurate extrapolation of D 0 , 
at least not to the accuracy of the present work. 

B. The Red and Violet Band Systems 

It is appropriate to address first the question of what level of treatment is re- 
quired to accurately compute the separations (T t ) and transition dipole moments. 
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In Table 3 we summarize the variation in these quantities with MRCI reference 
selection threshold and basis set. The T t values and transition moments are rather 
insensitive to both the basis set improvement and the threshold used in the selection 
of reference occupations. As basis set extension and the reduction of the reference 
selection threshold have opposite effects on the transition moment, both improve- 
ments would have to be incorporated together to improve our computed results. 
However, the small changes that occur as the level of treatment is improved do 
not warrant the expense of an MRCI(O.Ol) calculation in a [5s 4 p 3d 2/ lg] basis 
set. We therefore use an MRCI(0.05) treatment in the [4s 3p 2d 1/] basis for the 
determination of transition moments over the full range of r, and these results are 
given in Table 4. The computed moments are us-d together with the Klein-Dunham 
potentials to compute the Einstein coefficients and lifetimes. The experimental T t 
values were used in determining the transition energies. The variation in our MRCI 
A - X transition moment with r is in good agreement with the CASSCF results 
reported by Larsson, Siegbahn and Agren (1983). In addition, the r variation in the 
B - X transition moment is in good agreement with that inferred from experiment 
(Danylewych and Nicholls 1978; Arnold and Nicholls 1973). 


There is genera! agreement between previous experimental and theoretical de- 
terminations of the radiative lifetime of the B 2 E+ state. In Table 5, we give the 
Einstein coefficients for all transitions in the violet system ( B - X) for v' <9, where 


the intensity is at least 0.1% of the intensity of the 0-0 band. In Table 6, we com- 
pare our calculated lifetimes for v' <5 with previous determinations. The lifetime 
of the B state is determined principally by the B -> X decay: inclusion of the 
B _» A channel reduces the lifetime by about 1%. Our results show an increase in 
the lifetime with t>\ in agreement with the results of Cartwright and Hay (1982). 
The results of Larsson, Siegbahn and Agren (1983) show a decrease from »' = 0 
to v' = 1, but then an increase in lifetime with increasing v'. The high frequency 
deflection experiments of Duric, Erman and larsson (1978) are in general agree- 
ment with our calculation, showing an overall increase in lifetime with increasing v 
(within the reported error bars). Note that the value ascribed to v'=0 is the aver- 
age of measurements of 16 rotational levels ai d a measurement at the band head. 
The lifetimes of the individual rotational levels vary somewhat depending on the 
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degree of perturbation by the A 7 T\ state. Considering the increased uncertainty in 
the experimental lifetimes with vibrational excitation, the theoretical results should 
portray more accurately the variation in radiative lifetime with v The theoretical 
lifetimes are in overall good agreement with experiment (Duric, Erman and Lars- 
son 1978; Jackson 1974; Luk and Bersohn 1973) Based on the calibration studies 
in Table 1, we believe that our lifetimes would t e as much as 5% longer if the 0.01 
threshold selection were used. This would bring our lifetime for v'=0 into better 
agreement with experiment and with Larsson, Siegbahn and Agren (1983). How- 
ever, there appears to be very little remaining uncertainty in the description of the 

violet system. 

Einstein coefficients for v' <9 of the red system {A - X) are given in Table 7; 
again, only those bands with intensities of at least 0.1% of the 0-0 band are re- 
ported. These Einstein coefficients are used to compute the lifetimes of the A 3 II 
state which are summarized in Table 8, together with previous calculations and 
some experimental results. There is excellent agreement between our results and 
those computed by Cartwright and Hay (1982), while the lifetime results of Larsson, 
Siegbahn and Agren (1983) are shorter. Note that at the CASSCF level our A-X 
transition moment at 2.2 Oo is considerably larger than at the MRCI level— see 
Table 3. External correlation, therefore, increases the lifetimes. In fact, it is likely 
that our lifetimes in Table 8 are lower bounds, and that further refinements in the 
theoretical treatment, such as expansion of the CASSCF active Bpace and reference 
list to include 6 orbitals, would slightly (about 5%) increase the lifetimes. 

All of the theoretical treatments show a decrease in lifetime with increasing v\ 
and all predict the lifetime to decrease by about a factor of two between v =0 and 
r'=9. The only experimental results consistent with these predictions are those of 
Taherian and Slanger (1984), although their values show a somewhat more rapid 
variation of lifetime with t/. This may be contrasted with earlier experiments 
(Jeunehomme 1965; Katayama, Miller, and Bondybey 1979; Duric, Erman, and 
Larsson 1978), which show almost no variation in lifetime with v’. Since it is 
unlikely that the theoretical values would display the wrong variation with v', there 
would appear to be some systematic error in the earlier lifetime measurements. 
In fact, our study suggests that even the experimental lifetimes of Taherian and 
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Slanger (1984), which have the correct variation with v\ contain some systematic 
error or that there are other decay mechanisms perhaps arising from mtersystem 
crossing between the X and A states, that shorten the lifetime. 

It is disconcerting that all of the recent experiments that measure the lifetime 
directly (Conley et al. 1980; Nishi, Shinohara and Hanazaki 1982; Jackson, Payne, 
Halpern and Tang 1982; Taherian and Slanger 1984) obtain values in the range of 
6-9 /is for t/=0, if indeed our theoretical valui of 11.2 /xs is a lower bound. How- 
ever, Davis et al. (1986) and Arnold and Nicholls (1972) have determined oscillator 
strengths (/ numbers) for several of the transitions in the red system that support 
a longer lifetime. Their results are compared with our work, that of Sneden and 
Lambert (1982), and that of Taherian and Slanger (1984) in Table 9. The relative 
magnitudes of the experimentally derived / numbers agree well with theory. We 
have estimated the lifetime for t/=0 level of the A 2 II state from the / numbers of 
Davis et al. (1986) and Arnold and Nicholls (1972) by scaling our computed lifetime 
based on the ratio of the /oo values from experiment and theory. This results in 
lifetimes of 13.6 /xs and 12.1 /xs, which lie between our theoretical values and those 
derived from the solar data (Sneden and Lambert 1982). Although the lifetimes 
derived from solar data appear to be slightly too long, the variation in the lifetimes 
with vibrational quantum number are in excellent agreement with theory. This is 
consistent with the fact that the absolute lifetimes and f numbers derived from solar 
data are potentially less accurate than the relative values, as they depend somewhat 
on assumptions about the opacity of the sun. However, we feel that there is suffi- 
cient data supporting longer radiative Ufetim as for the CN red system to warrant 
further laser-induced fluorescence studies of the A 2 n state. 

IV. Conclusions 

The Do of CN is computed from CN-C( 3 P)+N( 4 S) and also from 
Cj}-_ f C“( 4 5)+N( 4 5) corrected using the experimental electron affinities of C and 
CN. The Do values computed using these two different approaches agree to within 
0.05 eV. These results are corrected for one- particle and n-particle space incom- 
pleteness by an extrapolation using equivalent calculations and experimental results 
for NO, C 2 and N a : our best estimate for D 0 is then 7.65±0.06 eV. Our determi- 



nation of the lifetime of the violet system is in agreement with previous results, 
from both theory and experiment. Our lifetime for the red system is more consis- 
tent with the interpretation of the solar data than most of the laboratory-based 
experiments. Further experiments are probab y required to fully understand this 
difference. Overall, our theoretical results for D 0 and the lifetimes of the red and 
violet systems almost certainly represent the best, in the sense of both accuracy 
and internal consistency, available set of data for CN. 
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Table 1. Spectroscopic constants for CN versus the level of correlation treatment 
using the [4s 3 p 2d 1/] basis set. ——————=======i 



’’e(A) 

Cc(eV) 

u> ( (cm *) 

CASSCF(7) 

1.176 

6.758 

2073 

CASSCF(9) 

1.182 

7.101 

2020 

SOCl(7) 

1.176 

7.160 

2034 

SOCI(7)+Q 

1.177 

7.154 

2028 

SOCI(9) 

1.182 

7.M6 

2016 

SOCI(9)+Q 

1.183 

7.1:69 

2011 

MRCl(0.0l)(9) b 

1.181 

7.:j21 

2020 

MRCI(0.0l)(9)+Q 

1.183 

7.: 106 

2006 

Expt. 

1.1718 


2068.6 


° The number of electrons correlated is given in parentheses. The constants are 
based on a quadratic fit in 1/r to the theoretical energies on a 0.1 Oo grid. 

6 The MRCI is based on 52 reference occupations at the 0.01 threshold. 
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Table 2. A comparison of the Do values of CN, CN 

, N 2 , C 2 and NO. 

Basis 

MRCI MRCI+Q 

MCPF 

CN (0.01 reference selection) 




[4s 3 p 2d If] 

7.32 

7.31 

6.96 

[5s 4p 3d 2/ Iff] 

7.53° ,b 

7.52 

7.18 

A 

0.22 

0.22 


[5s 5p 4 d 3/ Iff] 



7.19° ,c 

best Estimate* 1 

7.65±0.06 



CN - (no reference selection) 




[4s 4 p 2d 1/] 

0.83 

9.82 


[5s 5 p 3d 2 f Iff] 

10.05 b 

10.06 


A 

0.22 

0.21 

0.24 

N 2 (no reference selection) 




(4s 3 p 2d If] 

9.23 

9.24 

8.85 

[5s 4 p 3d 2/ Iff] 

9.58 fc 

9.60 

9.11 

A 

0.35 

0.36 

0.26 

Expt. 

9.76 



error in [5s 4p 3d 2/ Iff] 

0.18 

0.16 


Cj (0.01 reference selection) 




[4s 3p 2d 1/] 

5.94 

5.90 


[5s 4p 3d 2/ Iff] 

6.11 

6.07 


A 

0.17 

0.17 


Expt. 

6.17* 



error in [5s 4p 3d 2/ Iff] 

0.06 

0.10 


NO (no reference selection)-^ 




[4s 3p 2d 1/] 

6.16 

6.26 

5.65 

[5s 4p 3d 2/ Iff] 

6.39 

6.50 

6.08 

A 

0.22 

0.24 

0.23 

Expt. 

6.61 
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0.23 


0.12 


error in [55 Ap 3 d 2/ Iff] 


• Bond length not optimized (r=2.2 a 0 used), MCPF calculations in the 
[5s Ap 3 d 2/ Iff] basis set suggest this introduce! a 0.002 eV erroT in Do. 

b One o orbital was deleted in the calculation rear r e ; for N 2 at the MCPF level 

this reduced Do hy 0.008 eV. 

C xhree a orbitals were deleted in the calculation near r«. Based on corresponding 
results for N 2 , this should decrease Do by less t nan 0.03 eV. 
d See text. 

* Thermochemical value of 6.17±0.11 eV from Kordis and Gingerich (1973). 

/ Computed at experimental bond length: the experimental u> e is used in the cal- 
culation of Do- 
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Table 3. A comparison the separations and moments with level of treatments 



Separations(cm' *) 


[4s 3p 2d If] basis 

Method 

j> 

1 

X 

A(B - X) 

CASSCF 

12473 

26302 

MRCI(0.05) 

11064 

26317 

MRCI(0.05)+Q 

10960 

25594 

MRCI(0.025) 

11283 

26069 

MRCI(0.025)+Q 

11148 

25716 

MRCI(O.Ol) 

11222 

26026 

MRCI(0.01)+Q 

11155 

25891 

[5s 4p 3d 2 f ly] basis 

MRC1(0.05) 

11029 

26336 

MRCI(0.05)+Q 

10918 

25589 


Transition Moment (a.u.) 


Method Moment (A - X) Moment(B - A) Moment(B - X) 

[4s Zp 2d If] basis 


CASSCF 

0.33528 

0.09175 

0.74638 

MRC1(0.05) 

0.30076 

0.08087 

0.69961 

MRC1(0.025) 

0.29866 

0.08173 

0.68642 

MRCl(O.Ol) 

0.29688 

0.08165 

0.68353 

[5s 4p 3d 2/ \g] basis 

MRCI(0.05) 

0.30086 

0.07928 

0.70314 


° All calculations were performed at r— 2.2 ao, using reference lists determined for 
this r value. Hence the 0.05 threshold values differ very slightly from the values in 
Table 4, which were determined using a combined reference list. 
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Table 4. Electronic transition moments( 
states of CN. 

a.u.) between the X 2 E + , A 2 II and f? 2 E + 

R(«o) 

A-X° 


B-A° 

B-X b 

1.6 

0.28650 




1.8 

0.30608 


0.12678 

0.56554 

2.0 

0.30913 


0.09975 

0.66644 

2.2 

0.29974 


0.08323 

0.68640 

2.4 

0.28179 


0.08206 

0.64249 

2.6 

0.25820 


0.09772 

0.54118 

3.0 

0.20217 


0.12138 

0.28700 

3.4 

0.14570 


0.09429 

0.18454 

4.0 

0.07628 


0.02798 

0.15445 

•The matrix 

element reported is 

< 2 n| 

> 


fc The matrix 

element reported is 

< 2 E + |*| 2 £ + >. 
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Table 5. Einstein coefficients for the £ 2 £+ - X 2 C + violet system of CN. 


Band 0 Strength 6 Band 0 Strength 6 Banc 0 Strength 6 Band 0 Strength 6 

5- 8 3.8 

6- 3 13 


0-0 

1000.0 

0-1 

72.6 

0-2 

3.9 

1-0 

112.1 

1-1 

824.1 

1-2 

120.7 

CO 

1 

H 

9.9 

2-0 

1.4 

2-1 

192.8 

2-2 

688.8 

2-3 

150.9 

2-4 

16.5 

2-5 

1.1 

3-1 

2.6 

3-2 

249.1 


3-3 

585.0 

3-4 

168.3 

3-5 

23.1 

3-6 

1.9 

4-2 

3.0 

4-3 

286.5 

4-4 

505.7 

4-5 

176.4 

4-6 

29.1 

4-7 

2.8 

CO 

1 

2.5 

5-4 

309.0 

5-5 

445.5 

5-6 

177.8 

5-7 

34.3 


6-4 

1.4 

6-5 

319.4 

6-6 

400.0 

6-7 

174.5 

6-8 

38.5 

6-9 

4.9 

7-4 

2.4 

7-6 

319.9 

7-7 

366.0 

7-8 

167.8 

7-9 

41.8 

7-10 

5.9 


8-5 

3.8 

8-7 

312.0 

8-8 

341.0 

8-9 

158.9 

8-10 

44.0 

8-11 

6.9 

9-6 

5.3 

9-7 

1.5 

9-8 

296.8 

9-9 

322.9 

9-10 

148.3 

9-11 

45.5 

9-12 

7.8 


“Upper state vibrational level given first. 

6 Einstein coefficients are given with respect to the 0-0 band. To convert to absolute 
values in seconds -1 multiply by 1.482 xlO*. 
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Table 6. Radiative Lifetimes r(ns) of the CN B‘ E + state. 


V 

This work 0 

lsA 0 - 6 

CH°- C 

Duric rf Jackson* Luk-^ 

0 

62.69(62.36) 

66.80(66.6) 

62.29(62.05) 

66.2±0.8 65.6±1.0 60.8±2.0 

l 

63.25(62.83) 

66.63(66.4) 

62.88(62.58) 

66.3±0.8 

2 

64.19(63.68) 

67.10(66.7) 

63.67(63.31) 

64.3±2.0 

3 

65.52(64.89) 

68.22(67.8) 

64.84(64.37) 

65 .6±3.0 

4 

67.23(66.44) 

69.72(69.2) 

66.42(65.79) 

68.2±4.0 

5 

69.32(68.33) 

71.43(70.7) 

68.42(67.63) 

67.3±5.0 


“Values in parentheses include also the contribution from the B - A channel. 
6 CASSCF calculations of Larsson, Siegbahn anc Agren (1983). 

C CI results of Cartwright and Hay (1982). 

rf HFD experiments of Duric, Erman and Larsson (1978). The value for t>=0 is the 
average of measurements for 16 rotational levels and a measurement at the band 
Lead, whereas for v >0 the measurements are only at the band head. 

“Laser induced fluorescence results of Jackson (1974); average of 14 rotational levels. 
•^Luk and Bersohn (1973). 
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Table 7. Einstein coefficients for the A*Tl - ■X -2 E f red system of CN. 


Band 0 Strength" Band 0 Strength" Banc 0 Strength" Band 0 Strength" 


6-5 

6-6 


78.4 

138.8 


8-7 

8-8 


191.8 

27.7 


0-0 

1000.0 

o 

1 

142.0 

0-1 

331.9 

4-1 

777.0 

0-2 

34.1 

4-2 

898.5 

0-3 

1.1 

4-3 

38.2 

1-0 

1158.0 

4-4 

209.4 

1-1 

78.5 

4-5 

6.7 

1-2 

277.7 

4-6 

44.1 

1-3 

59.3 

4-7 

5.6 

1-4 

2.8 

5-0 

50.9 

tO 

1 

o 

742.3 

5-1 

424.9 

2-1 

787.9 

5-2 

1012.3 

2-2 

23.2 

5-3 

551.7 

2-3 

147.1 

5-4 

7.0 

2-4 

66.2 

5-5 

196.6 

2-5 

4.4 

5-6 

1.0 

3-0 

354.3 

5-7 

28.6 

3-1 

1051.1 

5-8 

5.2 

3-2 

290.7 

6-0 

17.0 

3-3 

142.2 

6-1 

194.0 

3-4 

50.5 

6-2 

738.3 

3-5 

58.7 

6-3 

989.1 

3-6 

5.4 

6-4 

234.7 


6-7 

13.0 

8-9 

36.2 

6-8 

15.7 

8-10 

2.1 

6-9 

4.3 

8-11 

2.3 

7-0 

5.4 

9-1 

10.6 

7-1 

78.8 

9-2 

89.9 

7-2 

420.5 

9-3 

390.5 

7-3 

959.9 

9-4 

887.5 

7-4 

774.8 

9-5 

919.1 

7-5 

50.0 

9-6 

234.4 

7-e 

155.5 

9-7 

37.1 

7-7 

74.8 

9-8 

181.3 

7-f 

27.4 

9-9 

4.2 

7-9 

7.0 

9-10 

37.6 

7-10 

3.3 

9-12 

1.5 

8—0 

1.7 



8-] 

29.7 



8-2 

204.8 



8-3 

675.9 



8—4 

1020.0 



8— > 

488.8 




“Upper state vibrational level given first. 

"Einstein coefficients are given with respect to the 0-0 band. To convert to absolute 
values in seconds -1 multiply by 65.47. 
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Tkble 8. Radiative Lifetimes r(fis) of the CN A*U state. 


V 

Theoretical 


SL“ 


Expt 




This work 

CH fc 

lsA c 


Conley** 

TS e 

Duric^ 

Jackson 0 

0 

11.2 

11.1 

8.1 

14.2 

7.3±1.5 

8.50±0.45 


7.2±1.0 

i 

9.7 

9.6 

7.0 

11.6 


8.02±0.55 



2 

8.6 

8.6 

6.3 

10.2 


6.67±0.60 

3.8±0.5 


3 

7.8 

7.8 

5.7 

9.0 


>.50±0.45 

4.0±0.5 


4 

7.2 

7.2 

5.3 

8.0 


4.70±0.23 

4.0±0.5 


5 

6.7 

6.7 

4.9 

7.4 


4 OA + 0 ® 5 

4.2±0.5 


6 

6.3 

6.3 

4.6 

6.8 



4.4±0.5 


7 

6.0 

6.0 

4.3 

6.3 



4.4±0.5 


8 

5.7 

5.7 

4.2 

5.9 



4.5±0.5 


9 

5.5 

5.5 

4.0 

5.5 



4.3±0.5 


a From analysis of solar spectrum by Sneden and Lambert (1982). 
fc CI results of Cartwright and Hay (1982). 



'CASSCF results of Larsson, Siegbahn and Agren (1983), using molecular orbitals 
optimized separately for both states. 

d Con ley ti al. (1980) — average of results for 2 il 1 / 2 and 2 n 3 / 2 measured by laser 
excitation. 

'Taberian and Slanger (1984)— extrapolation of time decays to *ero pressure fol- 
lowing photodissociation of C 2 N 2 - 

J D urk, Erman and Larsson (1978)- high frequency deflection measurements. 
p Jackson, Payne, Haipem, and Tang (1982)-lsser induced fluorescence. 
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Table 9. Oscillator strengths of the CN red system. 


Band 

This work 

Davis 0 

AN b 

SL C 

TS d 

o 

1 

o 

2.361 

1.94±0.06 

2.19 

1.9±0.2 

3.10±0.16 

1-0 

1.911 

1.50±0.03 

1.66 

1.5±0.2 

2.30±0.16 

o 

1 

cs 

0.908 

0.76±0.02 

0.77 

0.84±0.08 

1.17±0.10 

3-0 

0.335 

0.28±0.02 

0.28 

0.33±0.05 

0.48±0.04 

3-1 

1.350 

1.12±0.07 

1.15 

1.2±0.2 

1.92±0.16 

4-1 

0.770 

0.74±0.05 

0.64 

0.66±0.1 

1.18±0.06 


“Davis it al. (1986) from comparison of absorpt on measurements made in a high- 
temperature furnace between red and violet systems. 

* Arnold and Nicholls (1972) from shock tube studies. These data are based on a 


CN Do value of 7.89 eV. 

c Sneden and Lambert (1983) based upon analysts of solar spectrum. 
^Taherian and Slanger (1984) 
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